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We report on new types of two-component one-dimensional dark solitons (DSs) in a model of a 
dual-core waveguide with normal group-velocity dispersion and Kerr nonlinearity in both cores, the 
coupling between which is dispersive too. In the presence of the dispersive coupling, quiescent DSs 
supported by the zero-frequency background are always gray, being stable with the out-of-phase 
background, i.e., for opposite signs of the fields in the cores. On the contrary, the background 
with a nonzero frequency supports quiescent black solitons which may be stable for both out- and 
in-phase backgrounds, if the dispersive coupling is sufficiently strong. Only DSs supported by the 
out-of-phase background admit an extension to the case of nonzero phase mismatch between the 
cores. 


Dark solitons (DSs) are fundamental modes in me¬ 
dia where signs of nonlinear and dispersive terms in the 
nonlinear Schrodinger (NLS) equation, governing evolu¬ 
tion of excitations, are opposite. DSs were predicted 
in Ref. [l| in the context of the mean-field theory of 
Bose-Einstein condensates (BECs), and obtained in the 
framework of the inverse scattering method in Ref. 0. 
Experimentally, DSs were created in various physical sys¬ 
tems, including optical fibers SI , and BEC Q. 

DSs exist also in coupled NLS equations, including 
models of dual-core optical couplers Q. In the latter 
context, the inter-core coupling may be dispersive, like 
the cores themselves [7, Q , the effect that was observed 
experimentally [9j. The dispersive coupling introduces 
new physics, as it links the temporal structure of optical 
pulses with their energy distribution between the cores, 
similar to the coupling of the translational and spinor de¬ 
grees of freedom in spin-orbit-coupled (SOC) BECs ta¬ 
in particular, the interplay of the SOC with the cubic self¬ 
attraction of the condensate allows to suppress collapse 
and leads to formation of stable two-dimensional (2D) 
bright solitons in the free space 0- A similar mecha¬ 
nism produces stable families of spatiotemporal optical 
bright solitons in a dual-core planar w aveg uides with the 
Kerr self-focusing acting in each core [131 ■ SOC in self- 
repulsive BECs in optical lattices supports vortex and 
half-vortex solitons [l3|. Similarity of mathematical de¬ 
scription of evolution of SOC BEC and light propagation 
in waveguides with dispersive coupling may allow trans¬ 
fer of many interesting concepts from the field of optics 
to matter wave systems, and vice versa. 


DSs in dispersively coupled waveguides were consid¬ 
ered recently, mainly from the perspective of the design 
of switching devices EU [l5| . DSs were also studied in 
SOC-BEC models [HMT where, unlike in the optical 
setting, a trap potential is an inherent part of the phys¬ 
ically relevant model. In the presence of the trapping 
potential DSs bifurcate from the first excited state of the 


potential in the linear limit 00- The respective qui¬ 
escent DSs represent nonlinear modes with zero intensity 
at the center [0] (black solitons, BSs). 

In this Letter, we introduce an essentially new DS 
species in the dual-core waveguides, for which the pres¬ 
ence of the dispersive coupling is a necessary condition, 
i.e., they disappear or become completely unstable when 
the dispersive part of the coupling vanishes. We also for 
the first time illustrate nontrivial shape transformations 
of DSs due to considerable phase mismatch between the 
cores. 

We consider a system of coupled NLS equations for 
scaled field amplitudes qi. 2 in the two cores with the nor¬ 
mal intrinsic group-velocity dispersion and self-focusing 
Kerr nonlinearity: 
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Here z and r are the_propagation coordinate and reduced 
time, respectively [ 23 ], k > 0 is the linear coupling co¬ 
efficient, S is the coupling-dispersion strength, and (3 ac¬ 
counts for the possible phase-velocity mismatch between 
the cores [0 . In the most general form, we are interested 
in DSs satisfying boundary conditions (j = 1, 2) 

t hm^ q i = ex P {* i bz +ujt- (—l) J (f>/2 ± x] } , (2) 


where Qi, 2 are the background amplitudes in the cores, 
(j) is the phase difference between them, b and ui are the 
propagation constant and frequency, and x is a phase 
shift across the DS. 

In this work, we focus on the settings with the in-phase 
(cf> = 0) and out-of-phase {(j) = 7r) backgrounds. Accord¬ 
ingly, we define cr = exp (i(f>) = ±1. Then, the substi¬ 
tution of the ansatz © in Eqs. © connects the prop¬ 
agation constant and phase-velocity mismatch with the 
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FIG. 1: (Color online) Field profiles in the first and second (top and bottom rows) components for (a) out-of-phase GS, (b) 
in-phase GS, and (c) out-of-phase BS at b = 4, 5 = 1.5, /3 = 0, as well as for (d) out-of-phase GS and (e) out-of-phase 
“nearly-black” soliton with b = 2, S = 1, /3 = 3.4 (the latter one is a DS which would be black in the case /3 = 0). Background 
oscillations in the BSs are not shown, as they may be removed by transformation qi i2 —► 51,2 exp(±i5r) in Eqs. ifljl. see Eq. 
©• Here and in other figures k = 1. 


asymptotic amplitudes of the background at t —» ± 00 : nents. Such solutions have the form: 


, = (Q1 + Ql){<™ + Q1Q2) 3 _ (Qi~ QpjvK - Q1Q2) 

2Q1Q2 2Q1Q2 

(3) 

We start the analysis by producing exact DS solu¬ 
tions for /3 = 0 (no inter-core mismatch) and zero fre¬ 
quency, w = 0. Then, in-phase (<r = 1) and out-of-phase 
(<7 = —1) DSs, which, generally speaking, may move with 
nonzero velocity v, are obtained as q \ j2 = e lbz Wi t 2 , where 


w 1 = (JW 2 = [*(<7(5 — v) + a g tanh(a g (r — vz ))\, 
a 2 = b — ok — (08 — v ) 2 


(4) 


(subscript g stands for “gray”, see below). These DSs ex¬ 
ist above the propagation-constant cutoff corresponding 
to a 2 > 0, i.e., b > (oS — v) 2 + ok. A representative fea¬ 
ture of the DSs in the system with the dispersive coupling 
(<5 7 ^ 0) is that the quiescent solitons, with v = 0, unlike 
the conventional BSs, keep finite “grayness”, which may 
be characterized by the intensity of each component at 
the soliton’s center, r g = <5 2 , as it follows from Eq. (0. 
Therefore, the quiescent DSs with zero frequency ui are 
called gray solitons (GSs) below. Examples of complex 
profiles, uq i2 = wir,2r( t) + *UTi,2i(i"), of the in-phase and 
out-of-phase GSs are displayed in Figs. 0a) and (b), re¬ 
spectively. 

A second representative family of solutions of Eqs. m 
at f3 = 0 corresponds to solitons with nonzero internal 
frequency ui. Such solitons, satisfying boundary condi¬ 
tions 0, can be found in the form gi j2 = e lbz+l “ T w;i j 2 , 
where w = <5 and ui = —S respectively for the families 
with in-phase and out-of-phase first and second compo- 


w\ = 0 W 2 = a b tanh(abr), a 2 = b — ok + 8 2 /2. (5) 

with o defined above. The intensity distribution of such 
solitons resembles that in classical black states with zero 
intensity at r = 0 (hence we call this family BSs and use 
subscript “b” to distinguish it from the GS family), but 
at the same time they reside on the background wave 
with nonzero phase tilt oSt. An example of an out-of- 
phase BS is shown in Fig.QJc). 

BS and GS families are represented in Fig. [2] by the 
renormalized energy flow, defined as U = U\ + [/ 2 , with 

Ui t 2 = (\Qi, 2 ? - I«T,2(t)| 2 ) dr. For /? = 0, Eqs. 

0 and 0 yield [7 gi b = 4a gi b- The figure also demon¬ 
strates alternation of stable (black) and unstable (red) 
branches within each subfamily. The instability may 
stem from the modulational instability (MI) of the back¬ 
ground, and from the instability of the localized DS’s 
core. The MI is amenable to analytical investigation at 
j3 = 0. In particular, one can show that at u = 0 the 
out-of-phase background is always stable, while in-phase 
background is unstable. At u> 7 ^ 0 the out-of-phase back¬ 
ground is also stable, while its in-phase counterpart is 
unstable at <5 2 < k, and becomes stable for sufficiently 
strong coupling dispersion, S 2 > k. This result is ex¬ 
plained by the fact that substitution <ji j2 = gi i2 exp (*<5r), 
see Eq. ©, transforms Eqs. 0 into similar equations 
with coupling constant k = k — <5 2 . The sign flip of k at 
<5 2 > k is tantamount to replacing the in-phase configu¬ 
ration by an out-of-phase one, which leads to the stabi¬ 
lization. 
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FIG. 2: (Color online) The renormalized energy flow of out- 
of-phase [marked “cp = tt" in panel (a)] and in-phase [ u rp = 0” 
in (b)] BSs and GSs (subscripts “b” and “g”, respectively), 
vs the coupling-dispersion coefficient, <S. (c) The amplitude, 
ab,d, and grayness, rt>,d, of the out-of-phase BSs and GSs vs 
5. In all cases b = 4, /3 = 0. Black (red) curves denote stable 
(unstable) branches. Circles in (a) and (b) correspond to the 
GSs and BS displayed in Figs, ©a, c) and (b). 


We now turn to the effect of the phase-velocity mis¬ 
match, /?, on the existence and stability of two above- 
mentioned families of DSs. In this case, Eqs. © pre¬ 
dicts that background amplitudes for two components 
are different. For small /? <C 1, using 0 , one ob¬ 
tains Qj = {b — ct ) 1 / 2 — (—l) J 'e + 0 (e 2 ), where e = 
ft[b — an) 1 ' 2 /[2(2kct — b)}. Considering the simplest case 
of <5 = 0 (generalization for S > 0 is straightforward) 
one can rewrite Eqs. © in terms of the normalized fields 
Pj = W j/Qj U = 1>2) : 

(1/2)Pj,tt + k(Q3— j/Qj)(Pj ~P3—j) ~ Qj (\Pj | — 1 )Pj = 0. 

( 6 ) 

We look for solutions of these equations in the form of 
Pj = p — evj, where p is the BS solution of equation p TT — 
2(|p | 2 — 1 )p = 0 and v\ ^(r) — > 0 at r — >■ ±oo. One can 
show that vi = —V 2 and these functions are real, while 
v\ satisfies linear equation Vi jTT + 2 (2 + a — 3p 2 ) V\ = 
4p (1 — p 2 ). The asymptotic form of this equation at 
r —>■ oo is given by vi jTT = 0 for a = 1, and Vi iTT —4v\ = 0 
for a = — 1. This shows that for /3 > 0 there may be 
no exponentially decaying solutions, i.e., no extension to 
/3 ^ 0 is possible, for the in-phase solitons with a = 1, 
while out-of-phase solitons with a = —1 may exist. 

Systematic results for j3 ^ 0 were produced in a nu¬ 
merical form. In accordance with the above conclusion, 
DSs exist solely with the out-of phase background, see 


examples of profiles of GS and BS in Figs. |T] (d,e), re¬ 
spectively. At large /3, the DSs of both types feature 
a double-well structure of the first component (weakly 
pronounced in these examples), and the local intensity 
never vanishes (hence, the notation BS serves here only 
for stressing that this family resides on the background 
wave exp(—i5r) and that it becomes black at /3 = 0). 
Further increase of /? reveals new shapes of DSs, with 
the first component featuring an elevation (rather than 
depression) against the background. 

Properties of the DSs for /3 ^ 0 are summarized in 
Fig. [3] Note a non-monotonous dependence of the renor¬ 
malized energy flow on /3 for both GS and BS families. 
The existence of the GS family is limited to a finite do¬ 
main, /3 < (3 cr . The respective derivative dU/dfi diverges 
at /3 = P cr . The existence of this family is also restricted 
to finite values of the coupling dispersion, S. On the 
other hand, no such bounds were found for the branches 
originating from the BS at /? = 0. 



P 

FIG. 3: (Color online) The renormalized energy flow of the 
out-of-phase near-BSs (DSs, which become black at /? = 0) 
and GSs (a), and the energy-flow shares in the two compo¬ 
nents, si, 2 , (c), vs the inter-core mismatch /3 at b = 2, (5 = 1. 
The dashed and solid curves in (b) show, respectively, the 
amplitude (a) and grayness (r) of GSs vs /3. All families de¬ 
picted in the figure are stable. Circles in (a) correspond to 
the solitons in Fig. 0d),(e). 

In Fig. (3][b) we observe that /3 ^ 0 results in dif¬ 
ferences of the soliton amplitudes and grayness in the 
components, in addition to the difference in the back¬ 
ground amplitudes. The most pronounced difference is 
observed for the grayness, with a relatively weak ampli¬ 
tude mismatch, almost in the whole existence domain, 
except for a narrow region near (3 CI , where the amplitude 
of the first (second) component abruptly decays (grows). 
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Shares of the renormalized energy in the DS components, 
Sj = Uj/U, which are displayed in Fig. E^c), exhibit a 
non-monotonous behavior: the share of the first compo¬ 
nent in GS initially grows with /3, but near /3 cr it suddenly 
decreases and renormalized power of second components 
starts to dominate. A smoother non-monotonous depen¬ 
dence s- 1 , 2 ( 8 ) is also featured by the near-BS branches 
that do not have any cutoff in /?. 

Stability of all types of DSs described here was thor¬ 
oughly analyzed by means of computing the respective 
eigenvalues, using equations for infinitesimal perturba¬ 
tions linearized around the stationary DSs. Results are 
summarized in Fig. [4] For 8 = 0, the out-of-phase gray 
solitons (<r = —1) exist below the line of empty circles, 
given by S 1 2 3 4 5 6 7 8 9 10 11 12 = b—an, as obtained from a g = 0, see Eq. 0 
with v = 0. GSs are stable for relatively small b. A com¬ 
plex instability domain appears as b increases. For large 
propagation constant values GSs are unstable for small 5, 
while the increase of the coupling dispersion eventually 
stabilizes them. The dependence of stability domains for 
the out-of-phase GSs on 8 at fixed <5 = 1 is shown in 
Fig. Up,). For small 6, out-of-phase GSs are stable in 
their existence domain, which is located below the line 
of empty circles. For larger b , they are stable close to 
the upper boundary of their existence domain. In-phase 
GSs are always unstable and their existence domain is 
not shown here. 

Stability analysis for the in-phase BSs on the plane 
(b,S) is presented in Fig. QJc), where the stabilizing role 


of the dispersive coupling is evident. Note that back¬ 
ground of such solitons is modulationally unstable be¬ 
low line of open circles. Above this line there is a do¬ 
main of instability of localized soliton core, but at large 
5 values such solitons become completely stable. Finally, 
Fig. [Ud) shows rather complex stability and instability 
domains for the out-of-phase BSs on the plane of param¬ 
eters (b,S). Here, increasing dispersion of coupling also 
stabilizes solitons. The domains of stability for BSs on 
the (6, 8) plane are not shown, but at S = 1 such solitons 
are always stable at b < 4.75. 

In conclusion, we have reported new types of two- 
component DSs (dark solitons) in the model of two-core 
waveguide with normal dispersion in both cores and dis¬ 
persive coupling between them. The inter-core phase- 
velocity mismatch was included too. Due to the presence 
of the coupling dispersion, zero-frequency background al¬ 
ways supports GSs (gray solitons), which are stable only 
in the case of the out-of-phase background, with opposite 
signs of the fields in the two cores. On the other hand, 
the background with a nonzero frequency supports BSs 
(black solitons), which may be stable for the in-phase 
background too, provided that the coupling dispersion is 
strong enough. Solely the out-of-phase background ad¬ 
mits the extension of DSs to the system with the phase- 
velocity mismatch between the cores. 
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FIG. 4: (Color online) Existence and stability domains for 
out-of-phase GSs in the plane of (fo,<5) at /3 = 0 (a), and in 
the plane of ( b,5 ) at 5 = 1 (b); for the in-phase BSs in the 
plane of ( b,5 ) at /? = 0 (c), and for the out-of-phase BSs 
in the plane of ( b,5 ) at /3 = 0 (d). Stability and instability 
areas are marked by “s” and “u”, respectively. In (a) and (b), 
lines of open circles show the upper boundary of the existence 
domain. In (c) and (d) there is no upper existence boundary, 
while the horizontal line of open circles in (c) indicates the 
border where the background for the in-phase BSs becomes 
modulationally stable, S 2 = 1, see the text. In-phase GSs 
are supported by the modulationally unstable background, 
therefore they are not shown. 















